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Abstract. This paper is devoted to the study of statistical properties of the 
greatest common divisor and the least common multiple of random samples of 
positive integers. 



1. Introduction 

For any given integer n > 2, let us denote by x[ n \ X^\ . . . a sequence of 
independent random variables uniformly distributed in {1, 2, . . . , n} and defined in 
a certain given probability space endowed with a probability P. For a concrete 
realization we may take the unit interval with Lebesgue measure and Borel a- 



algebra as the probability space, and for j > 1, the variable Xj whose value at 



X [n) 

u G [0, 1] is 1 plus the j-th digit of the expansion in base n of Li. 

We are interested in studying the probability distributions of the random vari- 
ables 

gcd(x[ n \...,X^) and lcm(X[ n) ,...,XW) for r > 2, 

i.e., the gcd and the 1cm of random r-tuples of integers. 

Dirichlct's basic and classical result asserts that the proportion of coprime pairs 
of integers in {1,2,..., n}, 

-o#{(«, j) : 1 < i,j < n; gcd(i, j) = l}, 

tends to 1/C(2) = 6/tt 2 as n tends to oo, which, in the probabilistic setting intro- 
duced above, reads: 



lim P(gcd(Xr,xW) = l) = 



C(2) 

See, for instance, [14], Theorem 332. 

The limiting behavior of the whole distribution of the gcd of pairs follows imme- 
diately from the above: for each integer k > 1, 

i^P {gcd{ xf\x^)=k) = ^. 
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The probability distributions and moments of the gcd and the 1cm of pairs of 
integers arc described asymptotically in the following two known theorems: 

Theorem A. a) The mass function of the gcd of pairs of integers satisfies 

(1.1) P (g cd(X<»U<»>)= t ) = ^ +0 (I±i^>) forl < k <n. 
b) The moments of the gcd of pairs of integers are given by 

(1.2) E(gcd(X 1 ("Ur ) )) = Hn) + C + O(^) ; 

(1.3) E(gxd(xj"\^r) = ^^[^^-l]+OK- 2 ln(n)), for q > 2; 

The estimate (1.1) follows directly from the usual bounds in Dirichlet's result. 
Estimates (1.2) and (1.3) appear, for instance, in a paper of Cohen (take gin) = 1 
if n = 1, and g(n) = elsewhere, in the notation of Theorem in page 168 of [8]). 
The constant C is recorded explicitly there. See also [9]. 

Theorem B. a) For < t < 1, the distribution function of the 1cm of pairs of 
integers satisfies: 

(1.4) P(lan(*<»>, *<">) < ,„•) = 1 - ^ g ' ->'".7 '-»'» + 0,(!^). 
b) T7ie moments of the 1cm o/ pairs o/ integers satisfy 

(1.5) E(lcm(x{ n UW)») = c( ^ g + + 2 ] )2 n 2<? + 0(n 2 ^ 1 ln(n)), for g > 1; 

Observe that the bound of the error term in (1.4) depends on £. (Throughout 
the paper, the notation O t means that the constant in the O-bound depends only 
on t.) 

The estimate (1.4) is more involved than the corresponding result for the gcd; it 
is due to Diaconis and Erdds [9]. Notice that the denominator j 2 in formula (1.4) 
is missing in the statement of Theorem 1 in [9] . 

Result (1.5) can be traced back all the way to Cesaro (see [5], page 248). See 
also Theorem 2 in [9]. 

In this note we will prove a number of asymptotic results (Theorem A' and 
Theorems 1-3) concerning the probability distributions (mass distribution and mo- 
ments) of the gcd and the 1cm of r-tuples of integers, for r > 3. The case of gcd is 
quite direct, but not so the case of 1cm, as we see later. 

Theorem A can be readily extended to higher moments. It is worth recording it, 
as we shall use these estimates elsewhere (see [11]). 

Theorem A'. Let r > 3. 

a) For 1 < k < n, 

(1.6) P(gcd(X^,... > ^)) = fc ) = F l- y + 0(-l : 

b) Let q be a positive integer. 
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bl) If l<q<r- 2, 

(1.7) E(g«K*M . . . , x^y) = + 0p (*EM) . 

b2) // g = r - 1, 

(1.8) E( gcd (If, . . . , IW)'- 1 ) = ^ + O r (l). 

b3) Finally, for q > r, 

(1.9) E(gcd(xf ) ,...,X("))«) =5 ri9 ri«- r+1 +O r> ,(n«- r li 1 (n)), 
where the constant D r q is given by 

°- = (, + i)c(, + i) SC) ( - i)W;fa - r+t+i) - 

The estimate (1.6) is straightforward, and it appears, with no bound on the error 
term, in Cesaro ([6], page 293). See also, for instance, [7], [16] and [18]. For the 
sake of completeness, we prove Theorem A' in Section 3, particularly of (1.9). 

Observe that bl) implies in particular that the mean of the gcd of an r-tuple of 
integers in {1, 2 . . . , n} has a finite limit for r > 3: 

lim E(gcd(At\...,Al"))Hfcil, 

reflecting the fact that, on average, the gcd of an r-tuple is quite close to 1, for 
moderately large r. Notice also that the constant D2. q reduces to 

(,liK(,U (3 ""- t( « + 11) - 

as in (1.3). We should mention that the asymptotic estimates of moments of gcd 
above are valid also for non-integer q, but we limit ourselves to the integer case. 

The random behavior of the least common multiple of r-tuples is subtler for r > 3 
than for r = 2. The random variable 1cm, (X^, . . . , Xr ) could be normalized 
in several manners; for instance, in terms of n r , of its maximum possible value 
lcm(l, . . . , n), or in terms of the product of the Xj's: 

lcm(x[ n \...,X^) \cm(x[ n \...,X^ n) ) lcm(x[ n) ,...,X^) 

_ _ _ _ or — 

n r ' lcm(l,...,n) ' ■ ■ ■ X^ n) 

In all three alternatives, we obtain a random variable with values in [0, 1]. We will 
focus on the first alternative (but see Propositions 4.8 and 4.9 for the third one). 
Recall that Theorem B claims that the sequence of variables 

_ lcm(x{">,Af } ) n>l 

converges in distribution to a random variable C with values in [0, 1] whose com- 
plementary distribution function is given by 

P( r > t ) - 1 V i-^l-mQ-t)) 
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and whose moments are given by E(£ 9 ) = grcfe^nis ; <Z > 1- 

To state our results about the distribution function of the 1cm, we introduce the 
following notation: for r > 2 and s > 0, denote by 

(1.10) A r {s) := {(xi, ■ • ■ , x r ) : < X\, . . . , x r < 1 , xi • • • x r < s} 

the part of the unit positive r-cube where xi ■ ■ ■ x r < s. We write f2 r (s) for the 
volume of A r {s). Observe that O r (s) = 1 for s > 1. For s < 1, 

(1.11) ^( S) = s ^ ln(1/S)fc 



^ fc! 

fc=0 

(see Lemma 2.2). For each r > 2, write 

p 

The constant T r is the asymptotic proportion of r-tuples of integers that are 
pairwise coprime (see [4] and [19]). Clearly, lim r _>. 00 T r = 0; in fact, it does so very 
rapidly: lim^oo T^ r ln(r) = e r , where 7 is Euler's constant, see [15]. The first 
values are: T 2 = 1/C(2) « 0.60793, T 3 sa 0.28675, T 4 rj 0.11488, etc. 



Theorem 1. Letr>3. Then, for Q < t < 1, 
(1.13) liminf P(lcm(X 1 (n) , . . .,lW) < in r ) > 1 - -i_ 

n i/vi \\-L / / / I <r* \ 



1 - n^tf- 1 ) 



C(r) ^ J r 

(1.14) limsupP(lcm(^ n) ,...,X(")) < tn r ) <l-T r V ^|^ . 

As O r (s) = 1 for s > 1, the series in (1.13) and (1.14), for each < t < 1, are 
actually finite sums (the range extends to those j such that f~ x < l/t). Notice also 
that the right hand side of (1.13) is a distribution function, with value as t — >• 0, 
and value 1 as t — > 1. The right hand side of (1.14) is not a distribution function, as 
it takes the value 1 — T r C,(r) as i — >• (see Figure 1 for a depiction of the case r = 3). 

Setting r = 2 in Theorem 1, we recover (1.4), since T2 = 1/C(2) and ^(s) = 
S (l-M*)). 

For the moments of the 1cm of r-tuples, we prove: 
Theorem 2. Let r > 3. For eac/i integer q > 1. 

( , 15) h ,HM^,.»±e. 

ME C.[x{-l..,#)') >Tr (H, + iH) 

ri->oo n r 9 ((J + 1)' 

Again, for r = 2, since T 2 = 1/C(2), we recover (1.5) of Theorem B. 
Theorems 1 and 2 do tell us that 

P(lcm(X<"\ . . „XW) > trf) x £ LzMj^) , 
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and 

E(lcm(x[ n \...,X^) 9 ) xn r «, 

but, asymptotically, they only provide upper and lower estimates for the distribu- 
tion function and moments of the 1cm of r-tuples of integers in the limit n — > oo. We 
have not been able to handle the combinatorics that would lead to establish precise 
asymptotic estimates for general r; but we have studied in detail the case r = 3 
and obtained the precise result contained in Theorem 3. To state it, we need to in- 
troduce some more notation. We will denote by co(rn) the number of distinct prime 
factors of m, and for each r > 2 we will write T r (m) for the arithmetic function 
given by T r (l) = 1 and 

p\m 

The function T r is multiplicative, and T r (m) < 1, for m > 1. The case r = 3, 

i/p 



T 3 H=n) 



+ 2/p' 

p\m 

will be of special interest. Finally, we shall denote by J the Dirichlet series: 

z — ' m s 

m— 1 

Theorem 3. a) ForO<t<l, 

r(n) Y (n) v (n)\ ^ , 3\ 



lim P(lcm(X} n; ,^" J ,XW) <in 3 



(1.19) =1 - T 3 2^-3 2^ —2 (l-0 3 (tjm)), 

j=l J m=l 

b) For integer q > 1, 

/ E(lcm(X 1 (n) ,^ n) ,^" ) )?) C(2g + 3) , 

1.20 lim 1 ' 2 3 ' ' = T 3 7 q J J(q + 2 , 

Incidentally, observe that the case t = of (1.19) implies the identity 

^ T 3 (m)3"( m ) _ 1 



^ m2 T 3 C(3) 



(1.21) J(2) = 

m— 1 

(see Remark 5.2). Figure 1 compares the exact value (1.19) with the lower and 
upper bounds given in (1.13) and (1.14) for r = 3. 

Finally, we discuss certain waiting time questions concerning sequences of suc- 
cessive gcds and lcms 

zi=xi, z 2 = gcd(xi,x 2 ), Z3 = gcd(xi, 0:2,0:3 ) 
wi — xi, w 2 — lcm(xi,x 2 ), w 3 = lcm(xi,x 2 , x 3 ) 

of integers xi,X2, ■ ■ ■ drawn uniformly and independently from {1, . . . , n}. 

The sequence (zj) decreases almost surely to 1, while the sequence (wj) increases 
almost surely to the number lcm(l, . . . , n). Their respective expected waiting times 
are dealt with in Section 6. 
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The paper is organized as follows: we will analyze the properties of the gcd of 
r-tuples in Section 3. Section 4 contains the proofs of the results concerning the 
lcm of r-tuples. The particular case r = 3 for 1cm is studied in Section 5. Finally, 
Section 6 contains the analysis of those waiting times related to gcd and lcm. 

Acknowledgments. We thank Fernando Chamizo from the Universidad Auto- 
noma de Madrid for some fruitful discussions. 



2. Preliminaries and notation 

Throughout the paper, for a vector x = (xx, . . . , x r ) of positive integers, we 
abbreviate gcd(x) = gcd(a^i , . . . , x r ) and lcm(x) = lcm(a;i, . . . , x r ). Also if (3 = 
(/3i, . . . , fir) is a vector of positive real numbers we economize and write x < (3 
to mean 1 < Xj < /3j, for j = 1, . . . r. Even further, for a positive number z we 
simplify and write x < z to mean that 1 < Xj < z, for j = 1, . . . , r. 

We shall encounter a few times in what follows sums of the form 



/(gcd(x)), 

where r > 1 and / is some arithmetic function. They can be readily seen to be 



x < n 



(2-1) £/(gcd(x))=5>*/)(j) 

x<n j— 1 

Here, /i denotes the Mobius function and the symbol * stands for the Dirichlet 
convolution. The expression above is valid also for r = 1, with the conventional 
understanding that gcd(j') = j, for any integer j > 1. 

Equation (2.1) is sometimes referred as Cesaro's formula (see [5], [6]). Its simple 
proof follows: 

Proof of (2.1). Using the properties of the Mobius function, we have that, for any 
arithmetical function F, 

n n 

(2.2) Yl ^(x) = ]>>(fc) £ F(*) = Y / Kk) F{ ~ k y)- 

x<n, gcd(x) — 1 A:— 1 x<n,fc|x fc=l y^ n /& 
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In our case, 



£/(gcd(x))=£/(d) £ 1 = £ f(d) Yl 1 

x<n d— 1 x<n,gcd(x)— d d=l y<n/d,gcd(y) — 1 

= E =E L7J E □ 

We will use the following notation for the summatory function Sa of an arith- 
metic function a: 

Sa(x) = a(j) for each x > 0. 

j<x 

The following arithmetic function: 

(2.3) D r (m) = m r - (m - l) r , m>l 

will appear several times in the paper; its summatory function is given by 

(2.4) SD r (m) = m r . 

Observe that D r (m) < rm r ~ 1 . 

For each integer q > 1, we denote I q (n) = n q , n > 1; its summatory function 
satisfies 

m m 9+1 

(2.5) SI q (m) = IM = E k " = 7X7 + " (m?) • 

fe=l fe=i q 

For each integer q > 1, denote by ip q the arithmetic function given by <p q (n) = 
(/i * Iq){n) (the so-called q- Jordan totient function). 

Lemma 2.1. For integer q > 1, t/ie summatory function of (f q satisfies 

' + O g (n 9 ) i/g>2, 



(2.6) Sp,(n) 



(9 + l)C(9 + l) 

■ n 2 + 0(rt Inn) ifq=\. 



. 2C(2) 

Notice that 951 is just Euler's ^ function and that the case q = 1 of Lemma 2.1 is 
just Dirichlet's theorem. 

Proof. Apply the well-known expression (see [1], Theorem 3.10) for the summatory 
function of the Dirichlet convolution of two arithmetical functions a and /3, 

(2.7) 5(a*/9)(s) = 5>(7) S/3(x/j) = £/3(j) Sa(x/j). 

j<x j<x 

and equation (2.5). □ 

The following is an elementary calculus lemma which shall be useful in the dis- 
cussion of distributional properties of 1cm. 

Lemma 2.2. For each r > 1 and for all s > 0, denote by £l r (s) the volume of 
the r-dimensional set A r (s) := {(xi, . . . , x r ) ■ < < 1 , Xl • • • X r < s}. 

TTien Q r (s) = 1 /or s > 1 and 

_ , . ^ ln(l/ S y . 

fl r (a) = s 2 / ., fors<l. 

3=0 J " 
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Proof. Observe that f2i(s) = s and that Q r (s) = s + f Q r -x(s/x) dx for r > 2. □ 

A simple change of variables gives that if < (3\ , . . . , f3 r < 1, then 

(2.8) Vol{(xi, . . . ,x r ) : < xj < = l,...,r, xi • ■ ■ x r < s) = BCl r (s/B) , 

where B — 17^=1 Pj- Observe that if B < s, the statement is obvious. 

Later on, Section 6, while discussing waiting times, we shall need the stan- 
dard Euler's extension to real argument of the harmonic numbers given by H n — 
Ylj—i 1/j for integer n > 1; and the standard approximation of H n by ln(?i), which 
we record in the following: 

Lemma 2.3. The parametric integral defined for any real a > by 

/>oo 

H(a) = [l - (1 -e- l ) a ]dt 
Jo 

satisfies 

H(a) = ln(a) + 7 + o(^— ^ , as a —> 00 . 
Besides, of course, for any integer n > 1, H(n) = H n . 

3. Probability distribution of the gcd of t-tuples 

In the present section, we prove Theorem A', which describes the mass function 
and the moments of the random variable gcd(x[ n \ . . . , Xr), for r > 3. 

Cesaro's formula (2.1), with f = 8%, where 5x(n) = 1 if n = 1 and is elsewhere, 
yields 

(3.1) P(gcd(xf\...,X(")) = l) = 45>(j) -f- 

3=1 J 

which readily gives that there exists a constant C r > such that 

P ( gcd(AT 1 (n) , . . . , X^ n) ) = 1) - -5- < C r - for any integer n > 1. 

<,(r) n 

As gcd(a;i, . . . , x r ) — k means that k\x\, . . . , k\x r and gcA{x\/k, . . . , x r /k) = 1, we 
deduce that there exists a constant C r > such that 



P(gcd(xW *(«>)= A}- ^ 



< C r ' 



which is statement a) of Theorem A'. 

Next, we turn to moments. Statements bl) and b2) are immediate consequences 
of the estimate (1.6): just write 

-(«) v(nh?\ _ _LV"^! . n (- V". k '' 



fc=l 



For q < r — 2, the sum above tends to £(r — <?)/C( r ) as n -> 00 with an error bound 
O r (l/n r_9_1 ), which in the worst case (r = q + 2) is O r (l/n); while the right hand 
side O term is, again in the worst case, O r (ln(n)/n). 

If q = r- 1, we get that E(gcd(A 1 (n) , . . . , X^f- 1 ) = ln(n)/C(r) + O r (l). 
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b3) The argument above would not work for q > r, since the purported error term 
turns out to be of the same order as the main term. Using Cesaro's formula (2.1), 
the identities (2.4) and (2.7), and Lemma 2.1, one writes 



E(gcd(x[ n \...,X^) q ) = i5> 9 tf) L?f = ^ £ D r {3)Sip q {\n/j\) 

3 = 1 3 3=1 

= — ( X1T77 ZiT Z^O ~ — 2-^0 

(9 + i)C(g + i)^ < U J q r 

As D r (j) < r j 1 * -1 and (7 > r > 3, this last error term is at most O q ^(n q ~~ r ln(n)). 
Finally, using the definition (2.3) of D r (j), and recalling again that q > r, we get 
that 



V Dr(j) _ I" ~ (j ~ 1 ) 7 ' - V 1 V f r V 1 \k+l -T-/. 

Z^ ,'g+l Z^ »g+l Z^ ? -g+l Z^ I U j v ' J 

J'=l i=l i'=l J fe=l V 7 

= 1 E jpiU ^ £ (I) <- 1 >' + '«' -' + * + !> + Ms) 

fe=l v 7 j=l J fe=l v 7 

and this proves (1.9). 

4. Probability distribution of the lcm of t-tuples 
Observe that 

af-ar rii<i<fe gcd(a i} a i; a fc ) 



(4.1) lcm(oi,...,Or) = „ „ 

rii<j g cd K a j) lli<i<fe<; gcd(a l7 a,-, a k , aj) 

Notice that the lcm is the product of the numbers, lcm(eti, . . . , a r ) = a\ ■ ■ ■ a r , if 
and only if they are pairwise coprime, that is, gcd(ai,aj) = 1 for each i ^ j. For 
r = 2, of course, there is no difference between coprimality and pairwise coprimality. 

4.1. Pairwise coprimality and equidistribution. For a r-tuple of positive in- 
tegers x, we write x G PC if the r components of x are pairwise coprime (that is, 
gcd(xi,Xj) = 1 for i ^ j). The following result was obtained by Toth and also by 
Cai-Bach (see [19] and [4]): 

Lemma 4.1. For each r > 2, 

(4.2) lim P({x[ n \ . . ..iW} £ PC) = J] (l - -Y^h + —) := T r . 



P 



When r = 2, the constant T2 is 1/C(2) and we recover (1.6). 
Extending the argument of Cai and Bach, one could prove that the PC r-tuples 
are, in fact, equidistributed, in the following sense (see the details in [10]): 

Lemma 4.2. Fix r > 2. Then, for any function f 6 C([0, l] r ), 



lim 



— V /( — ,...,—) =T r / f(u l> ...,u r )dui---du r 

x<ri,x£PC - 
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Taking / = 1, we recover Lemma 4.1. Let us record now two applications of 
Lemma 4.2 which will be useful in forthcoming arguments. Fix < t < 1 and f3 £ 
W with < ft, . . . ,ft. < 1. Let B = f]T =1 ft and define 6 = [0, ft] x ■ ■ • x [0, ft]. 
Choose/ = 1-BnA r (t)i where A r (t) the region given in (1.10). Then, thanks to (2.8), 

(4.3) lim — #{x < n[3;x £ PC; x x ■■ ■ x r < tn r } = T r BSl r (t/B) , 

n— too n r 

where 51 r is the function given in (1.11). 

For q > 1, take f(ui, . . . , u r ) = u\ ■ ■ ■ ■ lg. Then we get 

(4.4) lim V x\---x q r =T r . 1 . B q+1 . 

x<ri/3 vy ; 

xSPC 

Notice that actually the functions considered to derive these two examples are not 
continuous, as demanded by Lemma 4.2, but a standard approximation argument 
yields the results: for instance, for (4.3), consider 



to let e ->• 0. 



The following extension of Lemma 4.2 will be useful in Section 5. Fix an r-tuple 
a = (<2i, . . . , a r ) £ PC. We will say that x £ PC a if the components Xj are pairwise 
coprime and, additionally, gcd(xi,ai) = • • • = gcd(x,-,a r ) = 1. We refer, again, the 
reader to [10]. 

Lemma 4.3. Fix r > 2 and consider an r-tuple a = (ai, . . . , a r ) £ PC. Then, for 
any function f £ C([0, l] r ), 

lim ~7 ^ f(~ )•••)—) =T r T r (T\a^\ / f(ui,...,u r )dui---du r , 

x<n,xePC, j—1 • J i IJ - L \ 

where the function T r is given by (1.17). 

Taking / = 1, we obtain that the proportion of PC a r-tuples in {1, ... , n} r tends 
to T r Y r (Y[ T j = i a j) as n —> oo. Again, two special cases of interest: 

(4.5) lim — #{x < n f3; x £ PC a ; x x ---x r < tn r } 
n— >oo n r 

r 

= T r T r ( J[ aj^j B n r {t/B) , for t < 1 fixed, 
i=i 

and for q > 1, 

(4.6) lim —— V xl---a£ = T r T r (T]a j ) - f tt- B q+1 . 



x<ii/3,x£PC a 3 = 1 



4.2. Distribution function of the 1cm of r-tuples. In this subsection we shall 
prove Theorem 1. Fix r > 3. 
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4.2.1. Proof of (1.13). For < t < 1 and real z > 1, we introduce the following 
counting Functions: 

(4.7) N r {t,z) = #{x < z : x x ■ ■ ■ x r < tz r } 

(4.8) G r (t, z) = #{x < z : x x ■ ■ ■ x r < tz r , gcd(x) = 1} 

(4.9) L r (t, z) = #{x < z : lcm(x) < tz r } 

Our objective is to estimate 

P(lcm(xj" ) ,...,^")) <tn r ) = asrwoc, 

and we shall use for that purpose convenient auxiliary estimates for N and G. 

The following elementary lemma estimates the number of lattice points in the 
positive r-cube [1, z] r such that x\ ■ ■ ■ x r < tz r , in terms of the volume of the region: 

Lemma 4.4. For < S < t < 1, 



lN r (t,z) = n r (t) + O s (-) 
z r \zJ 



as z — > oo. 



Observe that the bound dependes on S and not on t. The next lemma estimates 
the function G r (t, z): 

Lemma 4.5. a) For < S < t < 1, we have that, as z —> oo, 

(4.10) G r (t,z) = ^z r + O s (z r - 1 ). 
b) For t > 1, 

(4.11) G r (t,z) = -^z r + 0(z r - 1 ). 

Proof. The estimate (4.11) is the case k = 1 of (1.6), as for t > 1, G r (t, z) counts 
coprime r-tuples. For the case t < 1, partitioning accordingly as the value of the 
greatest common divisor, we can write 

N r (t, z) = #{x < z : X\"-x r < tz r , gcd(x) = d} 

d<z 

= E ^ ^ : W • ' ' Wr < ^/d) r , gcd(y) = 1} = E G ^ 

By Mobius inversion (see, for instance, Theorem 2.22 in [1]), we deduce that 

G r (t, z) = ^2 ^ d ) N r{t, z/d). 

d<z 

Now, for t > 5, using Lemma 4.4 and r > 3, we get that, 

G r (t, Z) = J2 Kd) (J Or(*) + Oi (^)) = + ^l^ 1 ) 

= O r (i)^- + O a (^- 1 ), 
CM 

as claimed. □ 
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Our objective is to estimate L r (t, z). Let us write, again by partitioning, 
L r (t, z)=^2 #{x < z : lcm(x) < tz r , gcd(x) = d} 



d<2 



= ]T #{y < z/d : lcm(y) < tdT' 1 Q\ gcd(y) = 1} 

d<z 

>Y.#{y^ z l d - Vi — Vr< td^ 1 (^y, gcd(y) = 1} = °r(t d r ~\z/d) , 

d<z d<z 

where we have used that if gcd(x) = d and if we write xi = dyi, for i = 1, . . . , r, then 
gcd(y) = 1 and lcm(y) = lcm(x)/d. We have also used that y% ■ ■ ■ y r > lcm(y). 
The above inequality would turn into an equality if r = 2. 
Therefore, 



P(lcm(x[ n \...,X^) < triT) = > l-^GScf-^n/d). 



d<n 

To get the lower bound (1.13), we estimate the latter sum. For that purpose, we 
split it into two sums, depending on whether the first argument of G r is less than 
or equal to 1, or not: 

J2G r (t<T-\z)= £ G r (tcr-\Z)+ £ G r (t(f~ x , ^) := I + II. 

d<n d<n,td r - 1 <l d<n,td r ~ 1 >l 

For td r ~ Y < 1, thanks to (4.10), we can write 

where we have used, in the O term, that tdT^ 1 > t since d > 1. So the term I can 
be written as 

V ; d<n,id r - 1 <1 

where we have used that r > 3. 

On the other hand, using now (4.11) and r > 3, 

- e ^©'^(EGn^ e 

d<ra, td r - X >1 rf<« d<n,td r - 1 >l 

Adding up the expressions for I and II, 

Vr f*rfr-l n\_ nT V 1 + ^ V (^(trf r - X )-l) , n 

d<n ^ y 1 d<n ^ y 1 d<n. td r ~ 1 <1 



-»'(>- <k E ^«) + o, ( ,,->), 

from which inequality (1.13) is proved. 
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4.2.2. Proof of (1.14). We now obtain a lower bound of the complementary prob- 
ability by restricting to those r-tuples whose pairwise gcd are all equal and then 
partitioning according to the value of that common gcd. Fix < t < 1. Observe 
that 

n 

#{x < n : lcm(x) > tn 7 '} > #{x < n : lcm(x) > tn r ', gcd(xi, Xj) = k for i ^ j} . 

fc=l 

Writing Xi = ky i: for each i = 1, . . . , r, we get that y € PC and that lcm(x) = 
k • J/i • • • Vri a- n d conversely. So we can write 

n 

— #{x < n : lcm(x) > tn r } > — £ #{y < f : y x ■ ■ ■ y r > ^(f ) r , y € PC} 



k=l 



For fixed t, the sum has a finite number of terms, and the fc th one tends to 
j-pT r (1 — Q r (tk r ~ 1 )), as n -> oo, 



according to (4.3). So the whole sum tends 

v i-n r (tfc'- 1 ) 

r ^ k r ' 

tk T - 1 <l 

and (1.14) is proved. 



asn-> oo, 



4.3. Moments of the 1cm of r-tuples. In this subsection, we shall prove Theo- 
rem 2, with an argument akin to that in Theorem 10 of [13]. 

4.3.1. Proof of (1.15). The following lemma is a direct application of (2.2): 

Lemma 4.6. Fix q > 1 and consider the summatory function SI q (n) = J2j< n j q - 
Then, for r > 2, 



(4.12) 



gcd(y)=l 



Applying this lemma and the trivial estimate lcm(y) < y\ ■ ■ ■ y r we get 

^Ticm(x)« = £ lcm ( x ) 9 = ^E d9 E lcm W 



d < n x < n . 
gcd(x)=d 



gcd(y)=l 



(4.13) < ^ dV(fc)fc r9 [5/ g (^) 

This means that 

1 



E(lcm(^,...,X(«y <-i J2 



kd<n 



and we will get (1.15) by estimating this last sum. From (2.5), we deduce that 

, r(g+l)-l 



Mi) 



1 / n 



(g + i) r v&y 



Ma) 
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Plugging this into (4.13), the sum of leading terms is given by 

n r q 1_ C(r( g + l)-g) ^(,-lRl) 

while the sum of the error terms is seen to be O q . r (ji rq Adding this up, we 
obtain (1.15). 

4.3.2. Proof of (1.16). Arguing as in the proof of (1.14), 

n rq+r ^ / lcm(x) > rq ^_ r ^ y ^ < lcHl(x) — ^ ] n rq + r ^ ' ^ 1 ^ r 

x<n fc— 1 x<n /c— 1 y<n/k 

gcd(xi,Xj)=k,i^j yGPC 

According to (4.4) (with /3j = 1/fc), for each summand we have that 

}™L E Vl'"Vr= ( q+ l)r k r{q+l)-q 

y<n/k V ; 

yePC 



From the bound 

k q v- q k« /nyi (ny I 



y<n/k 

yePC 



dominated convergence and the fact that J2T=i l/k r{ - q+r >- q < +oo, we deduce 
^L^Tft E lcm(x)« = ^ C (r(, + !)-«), 

fe=l x<n vy ' 

gcd(a;i,a;j)=fe,i^j 

and therefore, 

hminf i E( lcm(x}"\ . . . , X™) = liminf * £ lcm(x)* > ^ ^ ' + g) 

x<n 

This proves (1.16). 

4.4. The logarithm of the 1cm. The structure of the 1cm exhibited in (4.1) 
invites to consider its logarithm, as it may be written in terms of sums of logarithms 
of gcd's of different lengths: 

r 

(4.14) In (lcm(ai, . . . ,a r )) = ^ln(aj) - ^ hi(gcd(a;, a.,)) 

3=1 i<j 

+ ^ hi(gcd(ai,aj,a fc )) 

i<j<k 

Using this, we can prove the following: 
Proposition 4.7. For r > 2, 

(4.15) ^[E(Mlcm(A('^...,A("))))-^ln(^ = £ Q (-1)* ^ . 
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Proof. Observe that, by Cesaro's formula (2.1), 

E(ln(lcm(X^. ..,*(«)))) = Q E (ln(X^)) - QE(ln(gcd(X^,^))) 



= ^EMi)-E(I)(-D fc E(M*i n )o-)f " 1v '' 



3 = 1 k=2 v ' j=l 



r 
n 



3 = 1 k=2 x 7 j=l 

We have used the fact that (/i*ln)(j) = A(j), where A denotes the von Mangoldt's 
function (see Theorem 295 in [14]). For k > 2 fixed, 

J™ E A W ( LyJ I) k = E tf 1 = ~ ^ . 

where we have used the trivial estimate A(n) < ln(n) and the standard expression 
of the Dirichlet series of the function A (see Theorem 294 in [14]). □ 



By the way, from (4.15) and Jensen's inequality one gets that 

E(lcm(X 1 ( " ) ,...,^"V) im . . 
liminf -i i-J_J U. > e l(H(r)-r) 



where H(r) is the function defined in the right hand side of (4.15), a lower bound 
weaker than (1.16), but of the proper order. 

4.5. Alternative normalization of 1cm. As we have mentioned before, one may 
alternatively normalize 1cm by dividing it by the product of the numbers. With 
an argument similar, but simpler, to the one used to prove Theorem 1, one may 
derive: 



Proposition 4.8. For r > 2, and for every < t < 1, 
liminf P( J ' < t) > 1 

n— >oo 



, ■ - T - 

X l ■•'■X-r W J jr -i <l/t J 

lim .l UpP ( Z&) ZJfi ^V^ 1 -^ E ~r 



1 ■■■^■r j>--i < i/t J 
Actually, for r = 2, we have equality for every < t < 1 : 

, lcm(X 1 (n) ,X (n) 
lim P' Vi>. 



n— >oo 



V x^x™ - > C(2) .^ /t f 



Of course, the statement above for r = 2 gives the asymptotic behavior of the 
distribution function of 1/ gcd(x[ n \ X^). The limiting distribution is discrete: it 
assigns mass ^jjp to the point 1/k, for every integer k > 1; in contrast to the 
limit distribution in Theorem B, part (a), which has no point masses. 

For moments, we have: 
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Proposition 4.9. For r > 2, and integer q > I, 

W E (( ^■■■x^ ))- r -« r fo + 1 >-*)- 

Actually, for r = 2, we have equality for every integer q > 1 : 

lim E(( \ - 1 ' 2 , 7 ) )= — -C(g + 2). 

»^oo VI X (n) -J^ n) ' ' C( 2 ) 

5. The case r = 3 

The analysis of this case r = 3 rests on a particular partition of the triples of 
integers which is based in the following factorization lemma: 

Lemma 5.1. To any triple of integers [x, y, z) we may assign uniquely an integer D 
and triples of integers (a,b,c) £ PC and (u,v,w) € PG( c ,b,a) so that 

x = D (ab) u, y — D (ac) v, z = D (be) w. 

In fact, D = gcd{x,y, z), 

{a = gcd(x, y)/D, (u = x/(Dab), 

b = gcd(x, z)/D, and <v = y/(Dac), 
c = gcd(y,z)/D, [w = z/(Dbc). 



Moreover, 



lcm(x,y, z) — D(abc)(uvw) . 



The proof is direct; we just insist that (u,c),(v,b),(w,a) are required to be 
coprime couples. For pairs of integers (x, y) the analogous representation is x — 
Du, y — Dv, with D — gcd(cc, y) and u, v coprime. 

5.1. Proof of part a) of Theorem 3. Fix an integer n > 1 and < t < 1. We 

follow Lemma 5.1 and partition the required counting: 

#{1 ^x,y,z < n; lcm(x, y, z) > tn 3 } 

= £ H #\ v - n /( Dac ) ■ (u,v,w) e PC (CjM) , uvw>t j^\ 

D=x (o,6,c)6PC I w < n /(Dbc) 

Now, according to (4.5), the argument of this double sum, for fixed D and fixed 
a, b, c, is asymptotically 



a a T 3 T 3 (abc) J3 * 2 (l - tt 3 (tD 2 abc) 



' i ■ i \ - > i i n > t i 

\abc) 
and is bounded by 

1 



< n 3 



D 3 (abcy 
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Since 

E E ^p^ = c(3) E p^ = c(3)E^<+°°' 

D=l (a,b.c)ePC V y (a,7>,c)6PC V ' m=l 

dominated convergence gives the result. Here we have used that for any given 
integer m there are exactly 3"^ m - ) triples (a, 6, c) S PC such that m = abc. 

Remark 5.2. For t = 0, equation (1.19) reads T 3 C(3) J(2) = 1. This begs to be 
proved directly. Observe that the function Y 3(771) = T3(m)3 u '^ Tn ' is multiplicative, 
and that, for any prime p and any positive integer a, Ta(p a ) = Y 3 (p) = 3 ^^^ . 
We can write the Dirichlet series defined in (1.18) as a product over primes: 

T(\ V- ggH TT^i , T 8 (p) t 3 (p 2 ) u{ 3(p+l) \ 



m 15 

m=l 



SO 



^)=n(i+7^-rT)=n 



+ p + 1 

(p + 2)(p - 1)7 ~ 11 V (p + 2 )(p-l) 

The reader may check, from the definition of T3 in (4.2) and the Euler product 
expression for C(«) = ]1 P t=^t, that T 3 C(3) = rj p ^g' 

5.2. Proof of part b) of Theorem 3. Fix q > 1. Lemma 5.1 allow us to write 

oo 

]T lcm(x,y,zy=J2 E (^( a5c ) 9 EL. Q&c ( uww ) 9 )' 

x,y,z<n D = \ (a,bx.)ePC 

where for D and (a, 6, c) fixed, the corresponding sum Ylj-,, b extends over 
{(it, v, w) € PC, u < n/[D(a&)], u < n/[D(ac)], w < n/[D(bc)}}. 
By (4.6), each Y^D a b c * s seen *° asymptotically 

,9+1 

'(q+ l) 3 \Di(abc) 2 

and is bounded by 

Y' <n 3 "+ 3 ! 

^—^D; a,b,c 

Since 



Ek a ,, c -« 39+3 ^T 3 w ? -^(^ 



as n — > oo . 



1 



D 3 (abcy 



00 / 1 \ <j+l 00 Q«(fn) 

E E ^««( 75 3 W ) -C(2, + 3)E^<oo, 

D=l (a,fc : e)6PC v y m=l 

dominated convergence gives that 

E ^( X ,y, Z )^T 3 a 2 q + s)y 3 ^ +2 , 

x,y,z<n v y m=l 

as claimed in (1.20). 



18 



JOSE L. FERNANDEZ AND PABLO FERNANDEZ 



5.3. Case r = 3, with alternative normalization. We have for every < t < 1, 

iim p(Y- — ;~ J V; 1 ' ) / ) = r ; 



^1 ^2 ^3 D 2 m>l/t 

Also, for integer q > 1: 

lcm(xl n) ,4 n) ,X (n) 



lira E 



6. Waiting times 

Fix n > 2 and consider the following experiment: draw integers xi,X2, . . . uni- 
formly and independently from {1, . . . , n}, and calculate the sequences of successive 
gcds and lcms: 

zi=x 1 , z 2 = gcd(xi,x 2 ), z 3 = gcd(xi,x 2 ,x 3 ), . . . 

wi = xi, w 2 = lcm(xi, x 2 ), w 3 = lcm(xi, x 2 , x 3 ), . . . 

The sequence (Zj) is decreasing, and each Zj > 1, while the sequence (wj) is in- 
creasing and each Wj < lcm(l, . . . , n). 

For each of these random sequences we are interested in the random variable that 
registers the first time when they reach their respective limiting values. Again, the 
case of the 1cm is quite more involved than the case of the gcd. 

6.1. Waiting time for the gcd. For fixed n, consider the (decreasing) sequence 
(Zm^) of random variables given by 

2{ B) =l{ n) , zW=gcd(Z^ 1 ,X(«)) = gcd(A 1 (n) ,...,xW) form>2; 

Denote by T n the first time when the sequence (zi™') reaches the value 1. The 
variable T n takes values 1,2,... 

Lemma 6.1. a) For fixed n, P(Z™' = 1) — > 1 as m — > oo. 

b) The mass function of T n is given by 

n 

_ — > / 1 n \ m 
(6.1) P(T„ > m) = -Y,^ k \-[j:\) form>l. 

and, for each m > 1, 

1 



(6.2) lim P(T„ > m) = 1 - 

n-»oo C( TO J 

Notice that the case m = 1 of (6.1) reduces to P(T n > 1) = 1 — 1/n. The case 
m = 1 of (6.2) is then obvious. 

Proof a) Observe that P(z[ n) = 1) = 1/n. Recall from (3.1) that, for m > 2, 



n 

P(i#> = 1) = P(gcd(At\ ■■■,X^) = l) = —Y J tik) 



n" 

k=l 



so 



|1-P(ZH = 1) |<^J_<1. 

fe=2 

and therefore, P(z£ ] = 1) -> 1 asm->oo (for fixed n). 
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b) Observe that, for each m > 1, the events {T n < m} and {Zm^ — 1} coincide, 
and therefore, 



1 " 

P{T n > m} = 1 - P{7; < m} = 1 - — ^ M*) 



fc=i 



fc=2 



From here one deduces immediately that, for m > 2, 



P(T„>m)- 1- 



1 



C(m) 



= O 



ln(n) 



This gives (6.2) for each m > 2. 



□ 



The above result tells us that T n converges in distribution to T, where T is the 
random variable given by P(T < m) = \/C,(m). Formulas for the expected values 
of these variables, T n and T, are given in the following result: 



Theorem 6.2. 

OC 

(6.3) lim E{T n ) = 1-1] 



k=2 



fe- 1 



E(T) 



QC 

£(> 



C(m) 



The numerical value of E(T) is around 2,7052; by the way, the sequence E(T n ) 
is not increasing. 



Proof. Changing the summation order, 



E(r„) = £ p(t» > -) = 2 - i - [JJ) m 



m=Q 



m=2 fc=2 



2 - ; - E (^iJ) = 2 - ^ - E"(*)^ffij 

k=2 m=2 fc=2 nLfeJ 



Now observe that 



< 



and lim 



As X)fe>2 l/(fc(fc — 1)) — 1) by dominated convergence, 

where we have used Landau's classical result (see [17]) that X)fcLiA*(&)A = 0- ^ 
can be seen seen (again by dominated convergence) that this sum coincides with 



E(T) 



m=2 



C(m) 



□ 



6.2. Waiting for the 1cm. Consider now the (increasing) sequence of 
random variables given by 

W[ n) = x[ n \ wW = lcm(wi"! 1 ,xW)=lcm(xW,... ) XW) for m > 2; 

Let us denote again by T n the first time when the sequence (wl"') reaches its 
limiting value, lcm(l, ... ,n). 
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In the analysis of T n we will use the following notation: for each prime p < n, 
denote by "f P {n) the integer such that 

(6.4) p^ (n) < n and > n. 
In formula, J p (n) = [lnn/lnpj. Notice that 

lcm(l,...,n) = Y[ P lp(n) - 

Now define f3 p (n) as the positive integer satisfying 

(6.5) p^ {n) (3 p (n) < n and p^ {n) (f3 p (n) + 1) > n, 

so that /3 p (n) = [n/p' fp ^\ . Observe that 1 < (3 p (n) < p. 

For each prime p < n, we consider the set C p (n) of integers given by 

C P (n) = {p^ n \2p^ n \. . .,/3 p (n)p*>W}. 

It is immediate to check that, for fixed n, the sets C p (n) and C q {n) are disjoint if p 
and q are different primes. Finally, call 

C*(n) = {l,...,n}- \JC p (n). 

Observe that the event of interest, { lcm(x[ n \ . . . , li"') = lcm(l, . . . , ri)}, can 
be written as 

{for each p < n, at least one among x[ n \ . . . , belongs to C p (n)}. 

This observation allows us to rewrite the waiting time question as a "weighted" 
coupon collector problem: we draw coupons (independently) from {1, . . . and 
to complete the collection of interest means to get, at least, one coupon from each 
of the classes C p (n), p < n. There are 7r(n) different classes, where 7r(n) denotes the 
number of primes < n, each one of them having "weight" f3 p (n)/n. The coupons 
from the set C*(n) are useless for our objective. See, for instance, [12], [3], or the 
survey [2] for information about a variety of coupon collector problems. 

In this language, the variable T n registers the time when the coupon collection 
is completed. Observe that P(T n > I) = 1 if I < 7r(n). In general, 

{% >1}=\J Ml) . 

where A p (l) is the event in which, among the first I coupons drawn, none of them 
belongs to C p (n). Applying the inclusion/exclusion principle, we can write 

P(T„ >i) = Y^ P(A p (0) - J2 P(^(0ni,(l)) + - 

p<n p<Cq<.n 



E 

p<n P <q<n 



From this expression for the distribution function of T n , and following [12], we 
get a compact formula for the expected waiting time: T n '- 



Lemma 6.3. For fixed n, 



(6.6) E(T n )=n l-H^ 1 -' 



— t P p {n)\ 

p < n 



dt. 
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Proof. Write 



E(r„) = ^p(T„>o = E(E( 1 

1=0 p<n 

f3 p (n 



(3 p (n) 



1=0 



p^p 



(3 p (n) + f3 q (nW 



p<n 2=0 

n 



E EO 

p<q<n l—O 

l 



The identity (6.6) now follows. □ 

If we were to care just for the specific coupons {p lp ^} p < n (the first coupon in 
each class C p (n)), then the time T n to collect all of them will be, of course, longer 
than the time T n ] and, in particular, on average : E(T n ) > E(7;) (see later (6.16) 
for a precise comparison). In this case there are exactly n(n) coupons of interest 
(with weight 1) out of a total on n coupons and, therefore, see Lemma 2.3, 



E(T„) = 



l-na-o 



dt = n 1 - (1 - e"*)^' dt = nH n 



(n) 



The asymptotic size of E{7" n ) may be obtained by appealing to the elementary 
Lemma 2.3 and to the standard error bound on the Prime Number Theorem: 



(6.7) 



7r(n) 



ln(n) 



< C 



ln(ra) 



(6.8) 



valid for each n > 1 (C is an absolute constant). We may write: 

E(f n ) = nH n(n) = n(ln(7r(n)) + 7 + 0{\n{n)/n)) 

= n ln(n) — nlnln(n) + 717 + 0(n/ln(n)) . 

We would like to obtain an asymptotic expression like (6.8), but for T n - For that 
purpose, we introduce the following frequency counting functions of the /3 p (n): 

(6.9) Wi(n)=#{p<n:j8 p (n)=j} ) j>l. 

A few properties of these uij are in order. 

Lemma 6.4. a) ujj(n) = 0, if j > ^Jn. 

b ) E?=i w i( n ) = n ( n )- 



(6.10) 



c) For j >1 such that j + 1 < yjn, 
n 



"n 

7T| - 



i + i 



Tr(y/n). 



T ) < Wj -(n)<7r(j)- 

d ) J2%i 3^j{ n ) = J2 P <n A>( n ) ~ n\n(2), as n ->• 00. Actually, 

(6.11) -V &,(n) = ln(2) + G>f-^- T N ) , as n->oo. 

n Vln n) / 



e) For each j > 1, 



3.12) 



Wj (n) 



lim 

™-^oo 7r(ra) j(j + 1) 
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Observe that part d) means that asymptotically the proportion of useless coupons 
out of the total of n coupons (the set C*) is 1 — hi(2), about 30%. 

Proof. Statement a) is equivalent to f3 p (n) < yfn. To verify this, observe that 
from (6.4) we deduce that "t P (n) = 1 if \fn < p < n. In this range, if f3 p (n) > ^/n, 
we would get that pj3 p {n) > n, a contradiction with (6.5). Whenever 7 P (n) = a > 2, 
we have that f3 p (n) < p < n 1 / 01 < i/n. 
Claim b) is immediate. 

c) Let us introduce J = J n = [y/n\ — 1 (so that J + 1 < y/n). For 1 < j < J, 
primes p < n such that 

n n 

< p < - 

3 + 1 3 

satisfy 7 p (n) = 1 and P p {n) = j; all other primes p < n not included among these J 
classes satisfy p < -j^ < y/n. Thus (6.10) follows. 
For d), we write 

oo J 

E^^) = E^) = EiH7)- 7r (7Ti)] + E AM 

j=l p<n j=l J J P<7TT 

=E<j)— (7^) + E««)=E<j) + o(^), 

3 = 1 J P<7TT 3 = l 

where we have used summation by parts and the bounds 

£ *m*7tt'(tti)- (555) md "(ttt)-°I 

p<n/(J+l) v ; 

Now, from the error bound (6.7), we obtain 

J J 1 n 



i j-[3Mn/j) Vln(r 

and, since, 

1 /"^™ 1 / 1 

V - . . = / : , . . dx + Ol — 

rjjln(ti/j) 7i xln(n/x) Mn(r. 



j]n(n/j) Ji x\a(n/x) Vln(n) 
- =l v v ' 

=ln(2) 

equation (6.11) follows. 

Finally, (6.12) of statement e) follows by dividing (6.10) by ir(n) and invoking 
the Prime Number Theorem. □ 

The following lemma furnishes some precise asymptotic estimate for the fre- 
quency LU\. 

Lemma 6.5. 

WlW= 21^ + °(inV)J 
Proof. The error bound (6.7) readily gives that 

/ % ( n \ In- n \ 

7r(n) - 7T - = -:— r^r + O —5 ; 

v ; \2J 21n(ra) Vl n 2 (n)/ 

the result follows from (6.10) (for j = 1) and the bound Tr(y/n) = O ( ^ . □ 
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We are now ready to estimate E(7^). We start by rewritting formula (6.6) as: 

(6.13) E(T„) = n / l - TT (l _ e-«f<W)] dt. 

j<Vn 

By keeping just the factor corresponding to j = 1 in (6.13), we obtain the 
following lower bound: 

/>oo 

(6.14) E(T„) >n [1 - (1 - e -*)^W] = n^ l(n) , 

Jo 

where we have resorted to Lemma 2.3. 

For an upper bound: using that e - ^* < e~ 2t for j > 2, that J2 UJ j( n ) = ^i 71 ) and 
the identity 1 — xy = X — x + x(X — y), we may bound 



< n [!-(!- e-*) Ul W (1 - e- 8 )" 1 "'"" 1 '"'] dt 
Jo 

/•CO 

[l - (1 - p-*) wi ( n )] dt + n / (1 - e -*) Wl(n) [1 - (1 - e - 2t )^(»)-"i(")] ^ 

= n H Ul ( n) +n (1 - e-*)" l(n) [l - (1 - e - 2 *)*( n >-'"i-(»»)] dt. 
Jo 

Now, since 1 - x 5 < 5(1 - x), for < x < 1 and 5 > 0, we may further bound 
E(T„) <nH Ul(n) +n^(n)-u> x {n)) j (1 - e -*)"i(») eT 2 * dt 



(6-15) nH^ l(n) < E(%) < nH^n) t-/t prj- 



7r(n) — wi(n) r:r 7r(n) — Wi(n) 

= nH Ul(n) + n- — — < nH Ul{n) + n — 

(wi(ra) + + 2) wi(n)^ 

We have proved: 
Theorem 6.6. 77ie mean value of7~ n satisfies: 

7r(n) — wi(n) 
wi(n) 2 

Finally, 
Corollary 6.7. 

E(7^) = nln(ra) - nlnln(n) +n{^f- ln(2)) + 0(?i/hi(n)) . 
Proof. Using Lemma 2.3, the limit (6.12), for j = 1, and Lemma 6.5, we have that 

'ln(n) 



ff Wl(n) =Mwi(n)) + 7 + 0(^) 



Gb£o) +7 + °(ln^)) +0 (^) = ln Glni)) + ^ + °(h 1 R) 
Also, from the limit (6.12), for j = 1, we deduce 

7r(n) - Wi(ti) _ fhi(n) 



In 



wi(n) 2 

Combining these two estimates, we get the result. □ 

Observe that, as a consequence of Corollary 6.7 and (6.8), we deduce that 
(6.16) E(7T0-E(r) = nln(2) + o( I -^-y). 
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